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ABSTRACT 

We present astrometric measurements of eleven nearby ultracool brown dwarfs 
of spectral types Y and late-T, based on imaging observations from a variety of 
space-based and ground-based telescopes. These measurements have been used 
to estimate relative parallaxes and proper motions via maximum likelihood fitting 
of geometric model curves. To compensate for the modest statistical significance 

7) of our parallax measurements we have employed a novel Bayesian procedure 
for distance estimation which makes use of an a priori distribution of tangen- 
tial velocities, Kan, assumed similar to that implied by previous observations of 
T dwarfs. Our estimated distances are therefore somewhat dependent on that 
assumption. Nevertheless, the results have yielded distances for five of our eight 
Y dwarfs and all three T dwarfs. Estimated distances in all cases are > 3 pc. In 
addition, we have obtained significant estimates of V tan for two of the Y dwarfs; 
both are < 100 km s _1 , consistent with membership in the thin disk population. 
Comparison of absolute magnitudes with model predictions as a function of color 
shows that the Y dwarfs are significantly redder in J — if than predicted by a 
cloud-free model. 

Subject headings: stars: low-mass, brown dwarfs — astrometry 
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Introduction 



Determining accurate distances to brown dwarfs is important for a number of reasons. 
Firstly, distance is a vital quantity in establishing not only the space density of these 
objects, but also the luminosity function which can then be used to test models of star 
formation at the lowest masses. Secondly, distances allow the spectra of brown dwarfs to be 
placed on an absolute flux scale to provide more quantitative checks of atmospheric models. 
Thirdly, distances for the nearest objects allow us to construct a more complete view of 
our own Solar Neighborhood, allowing us to directly visualize the relative importance of 
brown dwarfs in the Galactic context. Sometimes, distance determinations produce results 
wholly unanticipated. For example, the J-band overluminosity of the T4.5 dwarf 2MASS 
J05591914— 1404488 (Figure 2 of Dahn et al. 2002) was unexpected despite its location near 
the J-band bump at the L/T transitio n (e.g., Looper et al. 2008), a feature thought to be 



associated with decreasing cloudiness ( jMarley et el. 



20101 ) . It has been s uggested, howeve r, 



that the overluminosity i s due to the presence of an unresolved binary ( IBurgasser et al 



2002 



Dupuy &: Liu 



20121 ). Similarly unexpected was the recent determination that young, 



field L dwarf s are often significan tly underluminous for their spectral types at near-infrared 



magnitudes (IFaherty et al. 



20121 ) . 



Some of the ear 



(120021) 



Tinney et al. 



i est pa r allax determinatio ns for brown dwarfs were by 



Dahn et al. 



( 120031 ) 



Wbaetal. 



( 120041 ) . once surveys such as the Two Micron 



All-Sky Survey (2MASS; Skrutskie et al. 2006), the Sloan Digital Sky Survey (SDSS; York 
et al. 2000), and the Deep Near-infrared Survey of the southern sky (DENIS; Epchtein 



et al. 1997) began to identi fy L and T dwa r 



programs by groups such as iMarocco et al. 



s in large numbers. More rece ntly, parallax 



(120101 ) and bupuv fc Liul (120121 ) have pushed 



astrometry measurements to the latest T spectral subclasses. With the discovery of Y 
dwarfs from WISE (Cushing et al. 2011; Kirkpatrick et al. 2012) we are now pushing these 



-4- 



measurements to even colder temperatures (jBeichman et al 



20121 ). In this paper we present 



distance and/or proper motion measurements for an additional eight Y dwarfs, along with 
three nearby late-T dwarfs, and present the first tangential velocity measurements for Y 
dwarfs. 



Observations 



Our set of objects includes all known Y dwarfs for which we have imaging data at a 
sufficient number of epochs for parallax an d proper motion estim ation. The exception is 



WISE 1828+2650, presented separately by 



Beichman et al. 



(120121 ). In addition, we have 



included three late T dwarfs from an i nvestigation of the low- mass end of the substellar 



mass function within 8 pc of the sun ( IKirkpatrick et al 



201l[ ). The complete sample is 



listed in the observing log shown in Table [TJ 

Each of these objects has be en observed a t two or three epochs by the Wide-field 
Infrared Survey Explorer (WISE; 



Wright et al. 



imaging observatio ns by the IRAC inst rument (IFazio et al. 



Telescope (Spitzer; 



Werner et al. 



12010 )) and at lea st four more epochs of 



20041 1 on the Spitzer 



(120041 ) ). the WFC3 instrument (IStraughn et al 



Space 



201lh 



of the Hubble Space Telescope (HST), and various ground-based observatories. The 
observatorie s and instruments used are list ed in the footnote of Table HJ and further details 
are given by 



Kirkpatrick et al. 



(12011 



20121) 



3. Astrometry Measurement Procedure 

Astrometric information was extracted from the observed images at the various epochs 
using the standard maximum likelihood technique in which a point spread function (PSF) 
is fit to each observed source profile. The technique was essentially the same as used 
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in 2MASS, the details of which are given by ICutri et al.l ( 120031 ). except that the source 



extraction results presented here were made using coadded images rather than focal-plane 
images. The positional uncertainties were estimated using an error model which includes 
the effects of instrumental and sky background noise and PSF uncertainty. The PSF and 
its associated uncertainty map were estimated for each image individually using a set of 
bright stars in the field, the median number for which was 14. Since the coadded images 
were Nyquist sampled or better, sine interpolation was appropriate during PSF estimation 
and subsequent profile fitting to the data. 

In order to minimize the systematic effects of focal-plane distortion and plate scale 
and rotation errors, our astrometry is based on relative positions using a reference star 
(or set of reference stars) in the vicinity of the object. For most objects we were able 
to find a reference star within ~ 10" common to all images except for those of WISE, 
due to the lower sensitivity of the latter. In order to incorporate the WISE data it has 
therefore been necessary to include bright reference stars which in general were much 
more widely separated from the brown dwarf (up to ~ 1 00"). Most of these were taken 



from the 2MASS Point Source Catalog (ICutri et al. 



20031 1 . In order not to let these stars 



significantly compromise the astrometry measurements from the more sensitive images with 
close reference stars, we used a hybrid scheme in which the bright stars were treated as 
secondary references, bootstrapped to the close reference stars using the images in which 
they were in common. 

The procedure is based on the following measurement model for the observed separation 
between the brown dwarf and reference star: 



ref 

a t - a it 



St - 5] 



ref 



, BD 



5 BD_ (r t + Ar t ) + z/ ,_^ 



, cat 



+ AO + v t - va 



(1) 
(2) 
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where at, St and ajf,5jf represent the extracted positions of the brown dwarf and ith 
reference star, respectively, estimated from the image at epoch t based on the nominal 
position calibration of that image; a^ at ,5 i cat represent the catalog position of the reference 
star, and Aaf 1 *, A£? at represent errors in the catalog position; u t , v' t represent the estimation 
errors for the brown dwarf, and Uu, v' it represent the estimation errors for the reference star. 
These estimation errors include the effects of random measurement noise on the source 
extraction as well as the residual effects of focal-plane distortion in the position differences. 
We assume that they can all be described by zero-mean Gaussian random processes. 

If we further assume that the Aaf*, Ac\ cat are described a priori by zero- mean 
Gaussian random processes with standard deviations substantially larger than the 
extraction uncertainties of the reference stars, then an optimal estimate of the brown dwarf 
position can be obtained from: 



at = a t + jj 2^ ( a i ~ a it > W 

$t BD = s t + ^ J2 ^ - ^ ef ) ( 4 ) 

where TZ(t) is the set of detected reference stars in the image at epoch t, and Nt is the 
number of stars in the set. 

The resulting estimates are included in Table [1] in the form of offsets from the nominal 
position of the brown dwarf at each epoch, and the set of reference stars used is given in 
Table [2j After having obtained d t BB and 5 t , the individual reference star catalog errors 
can then be estimated using: 



A~ cat I 1 / - BD , ref \ / K \ 

1 te£(i) 
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AV at = -r + ^£(& +^-^) ( 6 ) 

* te£(i) 

where £(i) is the set of all epochs for which the ith reference star is detected in the 
corresponding image, and iVj is the number of epochs in the set. 

These values can be applied as corrections to the catalog positions of the reference 
stars, enabling a corresponding time series of estimated brown dwarf positions to be 
obtained separately for each individual reference star via (CQ) and (j2J). The scatter in these 
estimates then provides a check on the assumptions regarding systematic effects such as 
focal-plane distortion and possible small proper motions of the reference stars themselves. 
We have included the effect of this scatter in the final quoted error bars in Table [TJ 



4. Estimation of Parallax and Proper Motion 



The measurement model incorporated the effects of parallax and linear proper motion, 
with approprate cor rection for the Earth-tra iling orbit in the case of Spitzer observations. 



The equations used (IKirkpatrick et al 



20111 ) were as follows: 



cos 5i(ati — cui) = Aa + fi a {ti — + iT tTig Ri ■ W (7) 
6t -St = AS + ns(ti - tx) - 7r tlig Ri ■ N (8) 

where ij is the observation time [yr] of the ith astrometric measurement, and Ri is the 
vector position of the observer relative to the Sun in celestial coordinates and astronomical 
units. N and W are unit vectors pointing North and West from the position of the source. 
Ri is the position of the Earth for 2MASS, SDSS, WISE, and HST observations; for Spitzer 
observations, Ri is the position of the spacecraft. The observed positional difference on the 
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left hand side is in arcsec, the parameters Aa and A5 are in arcsec, the proper motion fi a 
and fis are in arcsec/yr, and the parallax 7r trig is in arcsec. 

Maximum likelihood estimates, based on the assumption of Gaussian measurement 
noise, were made of five parameters: the RA and Dec position offsets of the source at a 
specified reference time, the RA and Dec rates of proper motion, and the parallax. The 
uncertainties w ere derived using the standard procedure for maximum likelihood estimation 



(IWhalen 



197ll ) using the positional uncertainties quoted in Table [3j The resulting estimates 
of proper motion and parallax and their associated uncertainties are given in Table [3J and 
the model fits with respect to the astrometry measurements are presented in Figures [U [21 
and[3j The chi squared values, x 2 > f° r the parameter fits in Table [3] are, for the most part, 
close to the number of degrees of freedom, A^f, indicating reasonably good modeling of 
position uncertainties. Formally, the probability of exceeding x 2 given N^i has a median 
value 0.29. 

The parallaxes that we present are, strictly speaking, relative parallaxes since no 
correction has been made for the small parallaxes and proper motions of the reference stars, 
most of which are relatively near jy. However, the expected correction for such effects is 
only ~ 2 mas (IDupuy fc Liull2012l ) which is at least an order of magnitude smaller than our 
typical astrometric uncertainties listed in Table El so in this error regime the distinction 
between relative and absolute parallaxes is unimportant. 

In order to check to what extent our parallax and proper motion estimates may have 
been affected by systematic effects of focal-plane distortion not properly modeled by the 
statistical assumptions of the previous section, we have compared the rms residuals of 
the above fits (obtained using multiple reference stars) with those obtained using a single 
reference star for each brown dwarf, and found that there was no significant difference. This 
suggests that whatever residual focal plane distortion errors exist, they are smaller than the 
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random errors of source extraction. 

We have converted our maximum likelihood estimates of parallax into most probable 
estimates of distance taking into account both the parallax measurements themselves 
and prior information. The latter includes an assumption that our objects are spatially 
distributed in a statistically uniform manner. Formally, that would imply that parallax 
values are distributed a priori as P{tt) oc 7r~ 4 ; the singularity at zero would then lead to 
difficulties in estimating the a posteriori most probable ir. Even though the zero parallax 



can be excluded on physical grounds, there is still a bias towards small va 
for S /N < 4, maximum likelihood parallax estimates become insignificant 



ues such that 



Lutz fc Kelker 



19731 ). Fortunately there is additional prior information to alleviate this problem; small 
parallaxes (i.e. large distances) can be excluded if they are inconsistent with the observed 
proper motion bas ed on an assumed velocity dispersion of the objects being studied 



(IThorstensen 



20031 ). 



With these considerations in mind, our estimates of distance, d, are based on the 
following assumptions: 

1. Our maximum likelihood parallax values, ttul, are distributed as Gaussians with 
standard deviation o^. 

2. Our objects are distributed spatially in a statistically uniform way, so that the a 
priori probability density distribution of d is proportional to d 2 . 

3. The distribution of tangential velocities of Y dwarfs in the solar neighborhood can 
be described by a Gaussian random process with mean and standard deviation V and 
oy respectively; we assume the values V — 30 km s " 1 and oy = 20 km s -1 respectively, 



representative of previous observations of T dwarfs (IFaherty et al. 



2009|). 



We then obtain the most probable distance, d, by maximizing the conditional 



-Im- 



probability density P(d\7iML, jUml), which from Bayes' rule can be expressed by: 

P{d\rr M L, A*ml) cx d 2 exp(-(/W - V) 2 /{2a 2 )) exp(-(7r M L - l/d)7( 2 0) (9) 

where /iml represents the magnitude of our maximum likelihood estimate of proper motion. 
Our distance estimates are presented in column 9 of Table [3J The error bars correspond to 
the 0.159 and 0.841 points of the cumulative distribution with respect to P(d\iiMh, /^ml)- 



5. Discussion 

As is evident from Table [TJ our observations represent a mixed bag in terms of 
telescopes (and hence spatial resolution) and time sampling since they were not specifically 
designed for astrometry, but rather for followup photometry of brown dwarfs detected by 
WISE. The quality of the observations was quite varied, and not always with sufficient pixel 
subsampling for the estimation of the high quality PSFs necessary for astrometry. In the 
case of Spitzer, for example, each observation consisted of a set of only five dithered images. 

In addition, the time sampling of the parallactic cadence is of key importance in 
the estimation of parallax. The ideal sampling involves observations at solar elongation 
angles of ±90°, and this is achieved by WISE, albeit with large position errors (typically 
~ 0.1 — 0.3"). These elongation angles are critical for an object on the ecliptic and less 
important at high ecliptic latitudes. For the parallax measurements described here, the 
worst example of poor sampling was WISE 1541-2250 for which all of the non-WISE 
observations were in one quadrant of solar elongation angle (see column 8 of Tabled]), so 
it is not surprising that the observations did not yield a significant parallax measurement. 
he previous measurem ent, corresponding to an estimated distance range of 2.2-4.1 pc 



(IKirkpatrick et al. 



20111 ). was based on even fewer observations and furthermore used 



position estimates for which the PSF errors were somewhat underestimated. Our present 
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result of > 6 pc therefore supercedes that estimate, but this object should clearly be 
revisited once a more optimal sampling of the parallactic ellipse has been obtained. By and 
large, however, there is a good correlation between the sampling cadence and the quality 
of the parallax estimate; future observations will be optimized both for image quality and 
cadence. 

Nevertheless, significant parallaxes {S/N > 3) have been obtained for five of the eight 
Y dwarfs and all three of the T dwarfs, thus providing distance estimates. Also, we have 
combined the latter with our proper motion estimates to yield tangential velocities, Kan- 
Of course, our estimated values, d and Vt an , are somewhat dependent on the assumed prior 
distribution of V ta , n in Eq. flU]), and the assumed similarity to the T dwarf distribution may 
not be valid if the Y dwarfs represent a significantly older population. In order to assess 
the sensitivity to this assumption, the distance estimates were repeated using a oy of 100 
km s _1 . It was found that for a parallax significance S/N > 4, the increase in oy led to no 
more than a 20% change (always in the positive direction) in d and hence Vt an . For lower 
values of S/N, Vt an becomes biased towards the a priori value, V, in Eq. (jHJ). Thus in 
Table [3] we quote Vt an values only for S/N > 4. Similarly, for S/N < 4 the reliability of our 
distance estimates is dependent on the validity of the assumptions regarding the a priori 
distribution of Vt an . 

On this basis we obtained significant values of Vt an for two of o ur Y dwarfs; both a re 



< 100 km s l , suggesting membership in the thin disk population (IDupuy &: Liu 



2012). 



Similar analy sis techniques, both in terms of the source extraction and parallax estimation, 



were used by 



Wright et al. 



(120121 ) to estimate the distance to the T8.5 object WISE 
1118+3125, inferred (with the aid of its observed common proper motion) to be a member 
of the £ UMa system, with excellent agreement with the known distance of that system. 

The distance estimates for the present sample, all of which are > 3 pc, have enabled 
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the estimation of absolute magnitudes. These indicate that luminosities plummet at the 



T/Y boundary ( IKirkpatrick et al 



20121 ) as illustrated by Figures S] and \5\ which represent 
updated versions of the absolute magnitude versus spectral type plots from the latter work. 
The steep decrease may at least partially account for the apparent scatter in absolute 
magnitudes of objects of the same spectral type, since in the YO regime an error of half a 
spectral type apparently corresponds to more than a magnitude difference in luminosity. 
More data will be required to make a definitive statement, however. 

The absolute magnitudes also provide valuable guidance for models in the ultra-cool 
regime. To this end we have compared our observational results with model-based and 
empirical predictions using plots of absolute magnitude as a function of color, as shown in 
Figure El The Mj verus J — H plot in the upper panel shows that the Y d warfs continue 



the tr end set by the L and T dwarfs based on the parallax observations of 



Dupuy fc Liu 



(120121 ) . A key feature is the turnover in the blueward progression of the color at Mj y 16 



at considerably redder J—H than predicted by cloud-free models (jSaumon &: Marley 



20081 ) 



as illustrated by the solid curve. Such beha vior is also apparent in the color-magnitude 



diagrams for cloud-free models presented by 



Leggett et al. 



fcoioh . The dotted/dashed 



curves in Figure [6] represent models incorporating the effect of clouds con taining various 



amounts of Cr, MnS, Na2S, ZnS, and KC1 condensates ( iMorley et al 



20121 ). as indicated by 



the sedimentation efficiency parameter, / se( j; lower values correspond to optically thicker 
clouds. It is apparent that these models can account at least partly for the relative redness 
of some of the J — H colors but they predict a blueward hook for temperatures below 400 
K which does not appear to be matched by the observations. Perhaps some of the scatter 
in J — H colors in Figure [6] might be explained in terms of a patchy cloud model; it is also 
possible that the inclusion of water clouds might improve consistency with the observations. 



Figure [6] does show reasonable consistency between observations and models based on 
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IRAC colors, i.e., M[ 3 . 6 ] and M[ 4 . 5 ] as a function of the [3.6]- [4.5] color. The only major 
discrepancy is that WISE 1828+2650, whose effective temperature is believed to be ~ 300 
K, falls at a location more indicative of 500 K on these plots. 

We thank C. Morley for providing the results of model calculations, and also the 
referee for very helpful comments. This publication makes use of data products from 
the Wide-field Infrared Survey Explorer, which is a joint project of the University of 
California, Los Angeles, and the Jet Propulsion Laboratory / California Institute of 
Technology, funded by the National Aeronautics and Space Administration. This work is 
based in part on observations made with the Spitzer Space Telescope, which is operated 
by the Jet Propulsion Laboratory, California Institute of Technology, under a contract 
with NASA. Support for this work was provided by NASA through an award issued to 
programs 70062 and 80109 by JPL/Caltech. This work is also based in part on observations 
made with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope 
Science Institute, which is operated by the Association of Universities for Research in 
Astronomy, Inc., under NASA contract NAS 5-26555. These observations are associated 
with program 12330, support for which was provided by NASA through a grant from the 
Space Telescope Science Institute. This paper also includes data gathered with the 6.5 
meter Magellan Telescopes located at Las Campanas Observatory, Chile. This research 
has made use of the NASA/IPAC Infrared Science Archive (IRSA), which is operated by 
the Jet Propulsion Laboratory, California Institute of Technology, under contract with the 
National Aeronautics and Space Administration, and also the SIMBAD database, operated 
at CDS, Strasbourg, France. 
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Table 1. Observing log and relative astrometry measurements 



Object Sp RA (nom) Dec (nom) Instr. Band Date Elong. Act cos S AS 



WISE J035000.32-565830.2 



Yl 57.501375 



-56.975006 



WISE J035934.06-540154.6 



Y0 59.892083 



-54.031703 



WISEP J041022. 71+150248.5 Y0 62.594667 15.046819 
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95 
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-0.743 
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-1 
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(0 
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78 


2 
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-1 
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(0 
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J 
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2 
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-0 
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(0 
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Spitzcr 
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2011-11-20 
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5 


-0.378 
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-1 


148 


(0 


131) 


Spitzcr 
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17 
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-1 
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(0 


075) 
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2010-01-13 


93 


4 
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-0 
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(0 
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WISE 
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-0 
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(0 
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-0 
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(0 
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9 


-0.468 
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-0 
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(0 
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W2 


2011-01-11 


95 


7 


-0.374 


(0.295) 


-1 


169 


(0 


301) 


Spitzcr 


[4.5] 


2011-01-19 


87 


5 


-0.500 


(0.196) 


-1 


284 


(0 


237) 


HST 


J 


2011-08-09 


-110 





-0.521 


(0.039) 


-1 


499 


(0 


040) 


Spitzcr 


[4.5] 


2011-11-20 


148 


8 


-0.564 


(0.194) 


-1 


969 


(0 


085) 


Spitzcr 


[4.5] 


2012-03-20 


26 


3 


-1.036 


(0.113) 


-2 


147 


(0 


178) 


WISE 


W2 


2010-02-16 


96 


3 


-0.001 


(0.188) 


-0 


052 


(0 


222) 


WISE 


W2 


2010-08-26 


-89 


2 


0.945 


(0.168) 


-1 


083 


(0 


193) 


WIRC 


J 


2010-08-29 


-92 


1 


0.997 


(0.172) 


-0 


965 


(0 


258) 


Spitzcr 


[4.5] 


2010-10-21 


-144 





1.411 


(0.066) 


-1 


538 


(0 


050) 


Spitzcr 


[4.5] 


2011-04-14 


39 


7 


1.175 


(0.134) 


-2 


621 


(0 


090) 



Table 1 — Continued 



Object Sp 


RA (nom) 
[°] 


Dec (nom) 
[°] 


Instr. 


Band 


Date 


Elong. 

[°] 


Aa 


cos <5 
'] 


AS 
["] 








Spitzcr 


[4.5] 


2011-11-19 


-172.8 


2.134 


(0.134) 


-3.894 (0.079) 








Spitzer 


[4.5] 


2012-03-20 


63.7 


2.320 


(0.126) 


-4.417 (0.132) 


WISE J053516.80-750024.9 >Y1 


83.820042 


-75.007019 


WISE 


W2 


2010-03-31 


-89.5 


-0.361 


(0.284) 


0.458 (0.317) 








WISE 


W2 


2010-09-28 


95.9 


0.172 


(0.266) 


0.693 (0.182) 








Spitzcr 


[4.5] 


2010-10-17 


77.1 


-0.006 


(0.151) 


1.221 (0.099) 








Spitzcr 


[4.5] 


2011-04-17 


-106.0 


-0.582 


(0.145) 


1.337 (0.160) 








HST 


J 


2011-09-27 


97.1 


-0.284 


(0.086) 


1.112 (0.036) 








Spitzcr 


[4.5] 


2011-11-20 


43.3 


-0.297 


(0.093) 


1.223 (0.137) 


WISEPC J140518.40+553421.5 YOp? 


211.326667 


55.572628 


WISE 


W2 


2010-06-08 


96.0 


-0.165 


(0.145) 


0.107 (0.198) 








WIRC 


J 


2010-07-26 


50.2 


-0.828 


(0.401) 


0.013 (0.216) 








WISE 


W2 


2010-12-14 


-88.8 


-1.388 


(0.161) 


-0.050 (0.284) 








Spitzcr 


[4.5] 


2011-01-22 


-128.5 


-1.829 


(0.130) 


-0.210 (0.125) 








HST 


J 


2011-03-14 


-180.0 


-1.723 


(0.118) 


0.155 (0.171) 








Spitzcr 


[4.5] 


2012-02-21 


-158.6 


-4.002 


(0.445) 


0.323 (0.203) 








Spitzer 


[4.5] 


2012-06-22 


82.1 


-4.862 


(0.148) 


0.260 (0.412) 


WISE J154151.65-225024.9 Y0.5 


235.465250 


-22.840358 


WISE 


W2 


2010-02-17 


-89.8 


0.206 


(0.532) 


-0.209 (0.708) 








WISE 


W2 


2010-08-15 


96.4 


0.041 


(0.175) 


-0.093 (0.173) 








FIRE 


J 


2011-03-27 


-127.5 


-0.959 


(0.196) 


-0.351 (0.198) 








Spitzcr 


[4.5] 


2011-04-13 


-144.3 


-1.019 


(0.120) 


-0.418 (0.140) 








NEWFIRM 


J 


2011-04-17 


-148.2 


-1.340 


(0.566) 


-0.430 (0.682) 








MMIRS 


J 


2011-05-14 


-174.4 


-1.260 


(0.113) 


-0.292 (0.128) 








Spitzer 


[4.5] 


2012-04-28 


-159.7 


-2.028 


(0.101) 


-0.568 (0.144) 


WISE J173835. 53+273259.0 Y0 


264.648083 


27.549758 


WISE 


W2 


2010-03-14 


-90.8 


-0.136 


(0.188) 


0.036 (0.204) 



Table 1 — Continued 



Object 


s P 


RA (nom) 

[°] 


Dec (nom) 

[°] 


Instr. 


Band 


Date 


Elong. 

[°] 


Aa cos 8 
["] 


AS 

["] 










WIRC 


J 


2010-07-26 


139.5 


-0.009 (0.074) 


-0.314 (0.060) 










WIRC 


J 


2010-08-29 


106.9 


-0.154 (0.271) 


0.046 (0.266) 










WISE 


W2 


2010-09-09 


96.3 


-0.064 (0.174) 


-0.309 (0.200) 










Spitzer 


[4.5] 


2010-09-18 


87.5 


-0.075 (0.103) 


-0.465 (0.110) 










HST 


J 


2011-05-12 


-148.5 


0.291 (0.056) 


-0.613 (0.048) 










Spitzer 


[4.5] 


2011-05-20 


-156.2 


0.336 (0.171) 


-0.544 (0.174) 










Spitzer 


[4.5] 


2011-11-26 


19.1 


0.276 (0.138) 


-0.982 (0.082) 










Spitzer 


[4.5] 


2012-05-12 


-149.2 


0.734 (0.124) 


-0.765 (0.109) 


WISEPC J205628. 90+145953.3 


Y0 


314.120417 


14.998147 


WISE 


W2 


2010-05-13 


-90.6 


0.027 (0.172) 


0.042 (0.167) 










WIRC 


J 


2010-08-29 


166.0 


-0.100 (0.168) 


0.316 (0.170) 










WISE 


W2 


2010-11-08 


96.1 


-0.015 (0.135) 


0.163 (0.144) 










Spitzer 


[4.5] 


2010-12-10 


63.8 


0.354 (0.148) 


0.437 (0.160) 










Spitzer 


[4.5] 


2011-07-06 


-142.0 


1.087 (0.229) 


0.853 (0.146) 










HST 


J 


2011-09-04 


160.5 


0.933 (0.035) 


0.810 (0.083) 










Spitzer 


[4.5] 


2012-01-06 


36.5 


1.152 (0.203) 


0.936 (0.231) 










Spitzer 


[4.5] 


2012-07-18 


-154.2 


1.889 (0.032) 


1.332 (0.042) 


WISEPA J025409.45+022359.1 


T8 


43.539375 


2.399750 


WISE 


W2 


2010-01-27 


94.9 


0.052 (0.085) 


-0.745 (0.119) 










WISE 


W2 


2010-08-05 


-90.7 


1.673 (0.139) 


-0.509 (0.100) 










WIRC 


H 


2010-08-29 


-113.7 


1.839 (0.144) 


-0.514 (0.181) 










WIRC 


J 


2010-08-29 


-113.7 


1.864 (0.120) 


-0.546 (0.189) 










Spitzer 


[4.5] 


2010-09-17 


-132.2 


2.164 (0.110) 


-0.536 (0.064) 










WISE 


W2 


2011-01-27 


95.1 


2.349 (0.193) 


-0.306 (0.231) 










Spitzer 


[4.5] 


2011-03-02 


60.8 


2.868 (0.068) 


-0.411 (0.126) 



Table 1 — Continued 



Object Sp RA (nom) Dec (nom) Instr. Band Date Elong. Act cos <5 A<5 



WISEPC J150649.97+702736.0 T6 226.708208 



70.460000 



WISEPA J174124.26+255319.5 T9 265.351083 25.888750 



Spitzcr 


[4.5] 


2012- 


03- 


07 


55 





5.504 


'0.049) 


-0.162 


(0.067) 


WISE 


W2 


2010- 


05- 


12 


95 


4 


-0.101 


(0.106) 


-0.010 


(0.118) 


WIRC 


H 


2010- 


08- 


29 


-9 





-0.479 


(0.083) 


0.064 


(0.065) 


WIRC 


,1 


2010- 


08- 


29 


-9 





-0.467 


(0.082) 


0.086 


(0.059) 


WISE 


W2 


2010- 


11 


18 


-89 





-0.488 


(0.119) 


0.069 


(0.220) 


Spitzer 


[4.5J 


2010- 


12- 


22 


-123 


5 


-0.573 


(0.069) 


-0.005 


(U.Udo ) 


Spitzer 


r a trl 

[4.5] 


2011- 


04- 


23 


114 





-0.707 


(0.077) 


0.990 


(0.136) 


Spitzcr 


r a trl 

[4.5] 


2012- 


01- 


23 


-155 


8 


-2.192 


(0.307) 


1.377 


(0.292) 


Snit^pr 


f4 51 


2012- 


05- 


25 


82 


4 


-2.278 


(0.148) 


2.130 


(0.115) 


2MASS 


.1 


2000- 


04- 


11 


-117 


8 


-0.069 


(0.138) 


0.188 


(0.104) 


SDSS 


z 


2004- 


09- 


16 


90 


1 


2.320 


;0.087) 


8.248 


(0.104) 


WISE 


W2 


2010- 


03- 


15 


-90 


7 


0.114 


'0.261) 


0.125 


(0.099) 


PAIRITEL 


H 


2010- 


04- 


09 


-115 


4 


-0.194 


(0.082) 


0.109 


(0.109) 


FanMt 


J 


2010- 


04- 


10 


-116 


4 


-0.136 


(0.038) 


0.231 


(0.058) 


FanMt 


H 


2010- 


04- 


10 


-116 


4 


-0.103 


(0.100) 


0.240 


(0.064) 


WISE 


W2 


2010- 


09- 


10 


96 


4 


-0.582 


(0.132) 


-0.596 


(0.184) 


Spitzer 


[4.5] 


2010- 


09- 


18 


88 


6 


-0.679 


(0.143) 


-0.334 


(0.179) 


Spitzcr 


[4.5] 


2011- 


05- 


20 


-155 


1 


-0.463 


(0.184) 


-1.317 


(0.200) 


Spitzcr 


[4.5] 


2011- 


11- 


20 


26 


3 


-1.265 


(0.052) 


-2.312 


(0.059) 


Spitzcr 


[4.5] 


2012- 


05- 


08 


-144 


2 


-1.132 


(0.088) 


-2.977 


(0.073) 



Note. — The columns represent the object name, spectral type from 



Cushing et al 



i20]J; 



Kirkpatrick et al 



I 2012), nominal RA and Dec position 



(J2000), the instrument (or telescope), band, UT date of observation, solar elongation angle, and the measured positional offsets (in RA and Dec) of 
the source from its nominal position. The key to the entries in th e Instrument column is as follows: 



WISE = Wide-field Infrared Survey Explorer (Wrig ht et al 



2010) 



HST = WFC3 camera on the Hubble Space Telescope ( Straughn et al 



Spitzer = Infrared Array Camera (IRAC) on Spitzer (Werner et al 



20lJ; 



2004); 



FanMt = Fan Mountain Near-infrared Camera (FanCam) (Kanncganti et al 



2009); 



FIRE = Folded-port Infrared Echellette at Las Campa nas Observatory i Simcoe et al 



MMIRS = MMT and Magellan Infrared Spectrograph l IMcLeod et al 



2008 



2010) 



2004); 



NEWFIRM = NOAO Extremely Wide-Field Infrared Imager at Cerro Tololo iS waters et al 



2009) 



PANIC = Persson's Auxiliary Nasmyth Infrared Camera at Las Campanas O bservatory I Martini et al 



PAIRITEL = Peters Automated I nfrared Imaging T elescope on Mt. Hopkins ( Bloom et al 



SDSS = Sloan Digital Sky Survey l lYork et al 



2MASS = 2 Micron All Sky Survey dCutri et al 



2004) 



2000) 



2006); 



2003) 



WIRC = Wide-field Infrared Camera on the 5-m Hale Telescope (Wi lson et al 



2003) 
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Table 2. Reference stars used 



Object Sp RA(rcf) Dec(ref) Separation Comment 







[°] 


[°] 


["] 




WISE 0350-5658 


Yl 


57.505458 


-56.976833 


10.4 








57.498042 


-56.985000 


36.6 


2MASS 






57.493500 


-56.961917 


49.6 


2MASS 






57.469708 


-56.975861 


62.2 


2MASS 






57.520292 


-56.993056 


74.8 


2MASS 






57.539625 


-56.972556 


75.6 


2MASS 


WISE 0359-5401 


Y0 


59.895458 


-54.033444 


9.5 








59.894458 


-54.021056 


38.7 


2MASS 






59.908333 


-54.042639 


52.3 


2MASS 






59.932417 


-54.040583 


91.1 


2MASS 


WISE 0410+1502 


Y0 


62.600125 


15.058056 


44.7 


2MASS 






62.580167 


15.039056 


57.6 


2MASS 






62.607333 


15.034722 


62.0 


2MASS 






62.618333 


15.049167 


82.7 


2MASS 






62.607083 


15.023306 


95.0 


2MASS 






62.622125 


15.043389 


96.3 


2MASS 


WISE 0535-7500 


>Y1 


83.824208 


-75.009278 


9.0 








83.793542 


-75.004417 


26.4 


2MASS 






83.811542 


-74.998583 


31.4 


2MASS 






83.771250 


-75.010028 


46.7 


2MASS 






83.769292 


-75.004389 


48.2 


2MASS 






83.823500 


-74.990444 


59.7 


2MASS 


WISE 1405+5534 


YOp? 


211.327083 


55.574778 


7.8 








211.343208 


55.584722 


55.0 








211.305583 


55.585333 


62.7 








211.273042 


55.574167 


109.3 


2MASS 






211.380417 


55.576639 


110.4 


2MASS 


WISE 1541-2250 


Y0.5 


235.464417 


-22.836833 


13.0 


2MASS 






235.466750 


-22.831694 


31.6 


2MASS 






235.473958 


-22.848833 


42.0 


2MASS 






235.477375 


-22.845306 


44.0 


2MASS 
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Table 2 — Continued 



Object Sp RA(rcf) Dec(ref) Separation Comment 









roi 


1 1 








235.482167 


-22.843861 


57.5 


2MASS 






235.467875 


-22.857306 


61.6 


2MASS 


WISE 1738+2732 


Y0 


264.643542 


27.547750 


16.2 


2MASS 






264.657750 


27.554444 


35.2 


2MASS 






264.640292 


27.535667 


56.5 


2MASS 






264.657750 


27.534028 


64.5 


2MASS 






264.640917 


27.530556 


72.8 


2MASS 






264.652417 


27.572083 


81.5 


2MASS 


WISE 2056+1459 


Y0 


314.117042 


15.000111 


13.7 


2MASS 






314.118667 


15.002556 


17.0 


2MASS 






314.120417 


15.007694 


34.3 


2MASS 






314.132917 


14.993361 


46.8 


2MASS 






314.106625 


14.999250 


48.1 


2MASS 






314.123833 


15.014750 


60.9 


2MASS 


WISE 0254+0223 


T8 


43.540792 


2.412833 


47.3 


2MASS 






43.537250 


2.386722 


47.5 


2MASS 






43.557958 


2.400333 


66.9 


2MASS 






43.512667 


2.395778 


97.1 


2MASS 


WISE 1506+7027 


T6 


226.736375 


70.461806 


34.5 


2MASS 






226.677125 


70.475250 


66.4 


2MASS 






226.750958 


70.443639 


78.3 


2MASS 






226.658042 


70.478833 


90.8 


2MASS 


WISE 1741+2553 


T9 


265.355375 


25.896583 


31.4 


2MASS 






265.341375 


25.893556 


35.9 


2MASS 






265.332083 


25.883611 


64.2 


2MASS 






265.360417 


25.905361 


67.1 


2MASS 






265.346958 


25.869194 


71.7 


2MASS 






265.339750 


25.905861 


71.7 


2MASS 
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Note. — Columns represent the object name, spectral type, the RA and Dec 
values of the associated reference stars, their separations from the object, and a 
comment column indicating which of the reference stars are in the Two Micron 
All-Sky Survey (2MASS) point source catalog. 
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Table 3. Parallax and Proper Motion Estimates 



Object Sp x 2 N d{ n a cos S /J,g tt Signif. d Vtan 

[«yr-i] ["yr- 1 ] ["] [sigmas] [pc] [km s" 1 ] 



WISE 0350-5658 


Yl 


14.22 


11 


-0.125±0.097 


-0.865±0.076 


0.291±0.050 


5.8 


q 7 +1.6 
d - '-0.4 


18 ±4 


WISE 0359-5405 


Y0 


13.02 


15 


-0.177±0.053 


-0.930±0.062 


0.145±0.039 


3.7 


r\ Q+l-3 

°- y -o.s 




WISE 0410+1502 


Y0 


11.53 


9 


0.974+0.079 


-2.144+0.072 


0.233±0.056 


4.2 


42 +i.2 


50 ± 10 


WISE 0535-7500 


>Y1 


11.80 


7 


-0.310±0.128 


0.159±0.092 


0.250±0.079 


3.2 


21 +13 




WISE 1405+5534 


YOp? 


9.16 


9 


-2.297±0.096 


0.212+0.137 


0.133±0.081 


1.6 


> 3.4 




WISE 1541-2250 


Y0.5 


15.21 


9 


-0.983±0.111 


-0.276±0.116 


-0.021+0.094 


< 1 


> 6.0 




WISE 1738+2732 


Y0 


15.22 


13 


0.348±0.071 


-0.354+0.055 


0.066±0.050 


1.3 


> 6.0 




WISE 2056+1459 


Y0 


6.64 


11 


0.881±0.057 


0.544+0.042 


0.144+0.044 


3.3 


7 5 +4.3 
'•°-1.8 




WISE 0254+0223 


T8 


5.67 


11 


2.578±0.042 


0.309±0.050 


0.185±0.042 


4.4 


4 9+1-0 


62 + 10 


WISE 1506+7027 


T6 


17.44 


11 


-1.241±0.085 


1.046±0.064 


0.310±0.042 


7.4 


3 4+0-7 
°-*-0.4 


27 ±4 


WISE 1741+2553 


T9 


9.90 


19 


-0.495±0.011 


-1.472+0.013 


0.176±0.026 


6.8 


r o+l.l 
°-°-0.6 


45 ±6 



Note. — Columns represent the object name, spectral type, chi squared of the parallax/proper motion fit to the estimated 
positions, number of degrees of freedom, proper motion in RA and Dec, the maximum likelihood estimate of parallax and its 
statistical significance, most probable distance (corrected for Lutz-Kclker bias), and the tangential velocity. Distance lower limits 
are based on a 2a criterion. Tangential velocities are quoted only for cases with parallax significance > 4, otherwise the Vtan 
estimate becomes strongly biased towards the assumed a priori mean value of 30 km s _1 . 
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Years since 1/1/2010 Years since 1/1/2010 




Fig. 1. — Proper motion and parallax fits to astrometry measurements of four of the Y 
dwarfs. Blue symbols represent observations from the ground and Low Earth Orbit (LEO), 
and red symbols represent Spitzer observations. The blue and red curves represent the 
corresponding model fits, respectively. The origins for the position offsets on the vertical 
(Motion) axes have been adjusted with respect to the values in Table HJ the AS and Aa cos 5 
values are relative to a constant position fit, so they are relative to the weighted mean of 
the a and 8. In addition, the AS values are offset for clarity by different amounts for the 
different plots. 
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Fig. 2. — Proper motion and parallax fits to astrometry measurements of the remaining four 
Y dwarfs. Color convention is the same as for Figure [D 
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Fig. 3. — Proper motion and parallax fits to astrometry measurements of the three T dwarfs. 
Color convention is the same as for Figure [U 
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T6 T8 YO 

Dwarf Spectral Type 



Fig. 4. — Absolute H magnit ude as a function o 



of the corresponding figure in 



Kirkpatrick et al 



spec tral type. This is a revised version 



( 120121 ) and include s the objects from th e 



present paper and the new parallax estimate f or WISE 1828+2650 ([Beichman et al. 



Kirkpatrick et al. 



20121 ). 



(120121 ) . which appears still 



The blue curve represents the relation used by 
to be an accurate representation of the absolute magnitude versus spectral type trend despite 
the fact that our results have been revised since the Kirkpatrick et al. paper was published. 
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Fig. 5. — Absolute W2 mag n itude as a function of spectral type. As with Figure H] it is 
taken from iKirkpatrick et al.l ( 120121 ) except for the inclusion of the objects from the present 
paper. It also includes WISE 1639-6847 (Tinney et al. 2012, submitted). 



31 




Fig. 6. — Absolute magnitude as a function of color. Large 



filled circles with error bars 



represent the objects from this paper, plus WISE 1 8 28+2 650 (IBeichman et al. 



included are the L and T dwarfs from 



20121 ). Also 



Dupuy fc Liul (120121 ). represented by open circles and 



small filled circles, respectively. For com parison, model cur v es are overplotted. The solid 



curve represents a cloud-free model from 



Saumon fc Marleyi (120081 ) . assuming g = 1000 m 



s 2 , K zz = 0. The numbers along this line represent the assumed values of effective tempera- 



Morley et al. 



(120 12 ) 



ture [K]. Also plotted (dashed/dotted lines) are four cloudy models from 
with the same assumed g and K zz , and with various values of the sedimentation efficiency 
parameter, / sed , as indicated. 



